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We provide a dynamical mechanism to generate localized features during inflation. The local 
' feature is due to a sharp waterfall phase transition which is coupled to the inflaton field. The 

key effect is the contributions of waterfall quantum fluctuations which induce a sharp peak on 
the curvature perturbation which can be as large as the background curvature perturbation from 
^SJ ' inflaton field. Due to non-Gaussian nature of waterfall quantum fluctuations a large spike non- 

» I \ Gaussianity is produced which is narrowly peaked at modes which leave the Hubble radius at the 

O |. time of phase transition. The large localized peaks in power spectrum and bispectrum can have 

interesting consequences on CMB anisotropies. 
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I. INTRODUCTION 



u 

4h' 

Simplest inflationary scenarios predict almost scale invariant, almost Gaussian and almost adiabatic perturbations 
and the observed cosmic microwave background (CMB) anisotropics are in very good agreements with these predic- 
tions [l[. However, with higher precision data expected to come from near future cosmological observations such as 
PLANCK, one hopes that some deviations from the predictions of single- field inflation may be observed. In particular 

■ a detection of non-Gaussianity would exclude a large class of inflationary scenarios. 

I — 1 ' There have been many attempts in the literature to generalize local features during inflation which may have non- 
trivial predictions on the power spectrum and bispectrum [34l3j. This is partly motivated from glitches in the CMB 
angular power spectrum at £ ~ 20 — 40. These local features may originate from models of high energy physics, 
, particle creation or field annihilation during inflation [3-113] • 

■ In many of these models local features are generated by temporal violations of slow-roll conditions during inflation. 
However, many of them are based on rather ad hoc mechanisms to violate slow-roll conditions. In this work, we provide 

C*") > a dynamical mechanism to generate local features in the power spectrum and bispectrum. The model consists of a 
t^J- | simple chaotic inflation potential, m 2 (j) 2 /2, coupled to a heavy field x- After (f> reaches a critical value <f> = ^> c (^> Mp), 
. x becomes tachyonic and a rapid waterfall phase transition takes place. In this sense the model is similar to hybrid 
CN| inflation [25|, |2g|. However, the key difference is that the potential is not vacuum dominated but the chaotic type 

inflaton potential is the main driving source of inflation. Furthermore, inflation continues even after the waterfall 

phase transition. 

In this model, inflation may be divided into three stages. At the first stage <fi > <fi c , inflation proceeds as in chaotic 
inflation. The second stage at which <j> < <j) c is fairly short, of the order of the Hubble time, and x becomes tachyonic 
and a waterfall phase transition takes place where x settles to its global minimum. This causes a small change in the 
inflaton effective mass, m 2 — > m\ = m 2 (l + C) where C << 1. Then the third stage proceeds as in chaotic inflation 
again. In order to bring the local feature into the observable windon of CMB, we demand that the phase transition 
takes place 50 — 60 e-folds before the end of inflation. 

The key point in our analysis is the evaluation of the waterfall quantum fluctuations and its contribution to the 
curvature perturbation on uniform density slices £ [27j . For studies of the waterfall contribution to £, in the context 
of hybrid inflation, see [28^49}. 

Due to the sharpness of the waterfall phase transition, the waterfall contribution in £ is narrowly localized around 
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the modes which leave the horizon at the time of waterfall transition. Furthermore, due to the intrinsically non- 
Gaussian nature of the waterfall contribution, which is in the form of £ ~ Sx 2 > a large spiky bispectrum is generated. 
We emphasize that this is a local dynamical effect intrinsic to the waterfall quantum fluctuations which are absent in 
other models, hence is an observationally testable feature genuinely specific to our model. 

The rest of paper is organized as follows. In Section [TT| we present our model and basic setup. In Section IIIII we 
calculate C, using the SN formalism. In Section [V] we calculate the dynamical and intrinsic non-Gaussianities in our 
scenario followed by conclusion and discussions in Section I VII Technical details about the power spectrum of the 
waterfall quantum fluctuations and higher order SN perturbations are deferred to appendices. 



II. THE MODEL 



In this section we present our setup in modeling local features during inflation. In our picture, we consider the 
simplest inflationary model, the chaotic potential m 2 (j> 2 /2, where we assume the mass of inflaton undergoes a sudden 
dynamical change at cf> — <j) c . We would like to see how a sudden small change in the inflaton mass can be modeled in a 
consistent dynamical way and then look for its observational consequences in CMB. One of the best known dynamical 
mechanism for inducing local feature during inflation is the idea of waterfall phase transition. For this purpose, we 
consider the following setup: 

m 2 ? A / o M 2 \ 2 g 2 2 9 



4 V A J 2 

where <fi is the inflaton field and x is the waterfall field. Formally potential (|2.ip is identical to hybrid inflation (2|| H|| . 
However, as we mentioned above, in our picture inflation is mainly driven by the field 4> so in our picture effectively 
inflation proceeds as in chaotic model with potential V ~ m 2 jj(j) 2 /2 with some effective mass m e /y which undergoes 
a small but abrupt change at (f> = <j> c . The waterfall field x is employed to induce this change in mass. 

In this model, like in chaotic inflation, inflation starts at <fi — (pi Mp, with Mp = (87rG) _1 for G being the 
Newton constant, so one can obtain 60 e-folds or more to solve the horizon and flatness problems. The waterfall field 
is very heavy during the first stage of inflation so it remains at its instantaneous minimum x = 0- Once the inflaton 
field reaches the critical value (f> — 4>c = M/g, the waterfall field x becomes tachyonic and quickly rolls down to its 
global minimum Xmin ~ M / A. The final stage of inflation after <f> > (j> c proceeds as in chaotic inflation again but 
with the effective mass of the inflaton, m+, given by 



where 



m\ = m 2 +g 2 ( X 2 ) = m 2 (1 + C) (2.2) 
q 2 M 2 

In our model, we assume the change in inflaton mass is small, C < 1. In order to bring this local feature into the 
observable scale, we assume that the short waterfall stage which lasts for about an e-fold or so begins at around 55 
e-folds before the end of inflation. 

It is instructive to look into different contributions to the potential before the phase transition where V — m 2 (f> 2 /2 + 
M 4 /4A. We have V(<f) c ) = m 2 <f) 2 /2(l + C/2). Having C <C 1 in our model corresponds to the assumption that the 
inflationary potential is dominated by the m 2 (j) 2 /2 term. This is in contrast to the standard hybrid inflation where 
the potential is vacuum dominated corresponding to C ^> 1. 

Before we proceed further we pause to compare our scenario to the one studied in [6j where the potential is still 
given by Eq. (|2.ip . However, in Q they are interested in the limit where C > 1 so inflation is mildly vacuum 
dominated. Hence their model is effectively a single field model where the waterfall quantum fluctuations can be 
neglected. In contrast, in the present paper we pay careful attention to the dynamics of the waterfall phase transition 
and calculate rigorously the curvature perturbations induced from the waterfall quantum fluctuations. For this to 
happen, we require 4> c > 10Mp. 

As usual the cosmological background is given by 

ds 2 = -dt 2 + a(t) 2 dx 2 , (2.4) 

where a(t) is the scale factor. It is more convenient to change the clock from the cosmic time t to the number of 
e-folds N via dN = Hdt where H = a/a is the Hubble expansion rate. Denoting the time when the waterfall phase 
transition takes place by N c , we further define n = N — N c . Hence n < for the period before the phase transition 
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and n > after the transition. We denote the end of the waterfall transition when \ has settled down to its global 
minimum by n — rif. With this notation, inflation in our model has the three stages: (a): n < 0, (b): < n < nt 
and (c): rif < n < N e — N c where N e is the time when the inflation ends. We set N e — N c ~ 55 so that the waterfall 
transition falls into the observable range. 

We are interested in the limit where the waterfall phase transition is fairly sharp, corresponding to n/ < 1 so 
the waterfall transition and symmetry breaking completion takes place in an e-fold or so. The key ingredient in our 
analysis is the dynamics of the waterfall quantum fluctuations during this short period and the curvature perturbations 
induced from them. To investigate this we will use the SN formalism [271 |50h52| . 

As mentioned before, we assume that the waterfall field is very heavy so it stays at x = and x = Xmin during the 
first and third stages. It is useful to introduce the dimensionless parameters a and f3 by 



m 
H 2 



6Mp 



6g 2 Mj 



M 2 
If 2 



QM 2 Mp 
m 2 4> 2 



6g 2 M 2 



(2.5) 



where in the last approximate equalities in both expressions it is assumed that <fi ~ <p c . The assumption that the 
slow-roll conditions hold during the first and third inflationary stages requires a <C 1. Demanding that ^ > 1 
for the waterfall field to be heavy, we require g 2 3> m 2 /Alp. On the other hand, from the COBE normalization 
we have m/Mp ~ I0~ 6 so g 2 ^S> 10~ 12 . Furthermore, we assume the onset of the phase transition and the time 
of sudden change in the inflaton mass to occur about 55 e-folds before the end of inflation so that it falls within 
the CMB observational window. In order for inflation to proceed long enough after the phase transition as in the 
standard chaotic inflation we require (j> c > 10Mp so g 2 < 10~ 2 M 2 /Mp. Combining this with g 2 ^> m 2 /Mp we get 
to 2 < ICT 2 M 2 or M > 10~ 5 M P . Finally, from the definition of C we conclude that g 2 /X < ICT 2 C. For example if 
we take C - 1CT 2 , we have g 2 /\ < ICT 4 and A = C- 1 g 2 M 2 /m 2 - 10 2 ff 2 M 2 /m 2 > ldV > ICT 8 . 



A. Inflaton dynamics 



As mentioned above during the first and second stage, inflation proceeds as in chaotic inflation with the potential, 



v-(4>) 



i 



M 4 

IT' 



(2.6) 



During the short second stage \ grows until the self-interaction term A% 4 becomes important. Then the self-interaction 
induces a large mass and x settles down to its local minimum. We denote this time by rif. After that, the waterfall 
field starts rolling across the valley determined by d x V((f>, \) — 0- This gives 



2 _ 2 
X Xmin 



X 



A 



(2.7) 



As seen from the above equation the local value of \ & t the third stage is dictated by the value of the inflaton field. 
Accordingly, during the third stage, the inflaton field experiences an effective potential given by 



1 



Solving the slow-roll equations of motion for (j), we obtain, for the first and second stages, 



4M P (N - N c ) = — 4M|, n = </>(n) 2 - 



I -Cln 



(2. 



(2.9) 



whereas for the third stage, 

8M 2 {N e -N) = -<fi + <f,(N) s 



1 + C 

C 



hi 



c 



1 + C0 2 



I - 



c p 

1 + C6 2 



C ^(N) - 4> A e 



0(C 2 



(2.10) 
(2.11) 
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As usual inflation ends at <f> — <fi e where the slow-roll conditions e, r\ <C 1 are terminated corresponding to 
\/2Mp. The slow-roll parameters are defined via 



e- p 



Ml fV 6 \ 2 , r0 V A , 



v = M 2 P ^t. (2.12) 
2 \V J p V 

In the present case, around the epoch of the phase transition, the slow-roll parameters are approximately expressed 
in terms of the parameter a as 

e~J7~|. (2.13) 

B. Waterfall Field Dynamics 

In this section we study the dynamics of the waterfall quantum fluctuations. One key point in our model, similar 
to hybrid inflation, is that before the waterfall transition x is very heavy so that it firmly stays at its local minimum 
X = and classically there is no background evolution of the waterfall field. Following the prescription in [29| we 
assume that for each horizon-size patch one would observe Sx 2 as a homogeneous classical background which varies 
smoothly over scales larger than the horizon scale. Thus on a given, sufficiently large scale, say the comoving scale of 
the present Hubble horizon size, one can calculate the mean value {5x 2 { n )) an d the fluctuation, 

A X 2 (n,x) = 5 X 2 (n,x) - (S X 2 (n)) , (2.14) 

where (Sx 2 ( n )} determines the homogeneous background while Ax 2 (to, x) gives rise to the curvature perturbations on 
super-horizon scales. 

With this understanding we now look into the dynamics of the background waterfall field and its quantum fluctu- 
ations in some details. The background waterfall dynamics is governed by 

x" + 3x'+(-(3 + 9 2 ^ + 3\jpj X = 0, (2.15) 

2 

where here and below the prime denotes a derivative with respect to to Neglecting the self-interaction term Aj^ 
during the second stage when the transition proceeds, an approximate solution for the background waterfall field is 
given by 29] 



X{n) ~ x(n = 0) exp 

where 



3 x 



(2.16) 



\otp- (2.17) 

Here it is worth noting the relations among the model parameters. Our model has four parameters: M, to, g and 
A. On the other hand, we have introduced six non-dimensional parameters: a, /3, e, to e x and C. Some of these are 
functions of time, but at leading order in the slow-roll approximation, we may consider them as constants. Then as 
a convenient set of independent parameters, we may choose e, C, e x and A. Then we have 

a e l 2 6Ae 2 C 

a = 3e, P = , r] = e , g = — 

2e 

M ~ — Ml , to ~ -. Ml . (2.18 

e e 

x x 

The dynamics of the waterfall field quantum fluctuations in Fourier space, neglecting the self-interaction term, is 
governed by 

5xZ + 35 X ^+(^-P + 9 2 t>) S Xk = 0, (2-19) 



2 H 2 r * H 2 
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where 

SXk = J 7^72 <5x(x)e~ jk - x = a wXk (n) + al k ^) . (2.20) 

Here ak and a k are the annihilation and creation operators, respectively, with respect to a suitably chosen vacuum 
and Xk(n) is the positive frequency function. Plugging the background value of (f> given by Eq. (|2.9p into Eq. (|2 . 19|) 
results in 

5x'i + 3<5x k + ( j| e ~ 2 " - 4 n ) = , (2.21) 

where k c is the comoving momentum of the mode which exits the horizon at the time of waterfall phase transition: 
k c = Ha(n = 0). As seen from Eq. (j2.2ip . e x H 2 measures the effective tachyonic mass of the waterfall field when the 
tachyonic instability develops. The assumption that the waterfall phase transition is sharp requires that e x 3> 1. The 
normalization of \k is determined by the canonical commutation relation, which gives 

6x^-^x^=0^. (2.22) 

For large and negative n, Eq. (|2.21j) can be solved by the WKB approximation by setting 6\k oc exp[So +S\ H ]. 

Taking account of the normalization condition (I2.22[) and choosing the standard Minkowski positive frequency in the 
limit n — > — oo, the result to first order in the WKB approximation is [32| 



H e -3n/2 

S Xk(n) = —7=^- 7TU ex P 



({k/k c ) 2 e- 2n ' - e 2 ^') ^ dri 



(2.23) 



After waterfall phase transition, n > 0, some (low k) quanta of the waterfall field become tachyonic and behave like 
classical random fields 31]. In particular there are modes which become tachyonic even before horizon crossing. In 
order to find the dynamics of the waterfall quantum fluctuations, it is convenient to divide the modes into large and 
small modes which we denote below by the subscripts L and S, respectively. Large modes are those which exit the 
horizon before the time of the phase transition n = and small modes are those which are sub-horizon at n = 0. 

The large modes cross the horizon sometime before the waterfall transition and after horizon crossing their profile 
decays like cx e -3 ™/ 2 as seen from Eq. (|2.23p . At the time of transition, n = 0, the WKB approximation fails, but 
apart from a factor of order unity, the amplitude at n = is approximately given by 

\5 X t(n = 0)\^-^=, (2.24) 

After the transition, n > 0, these modes follow the same equation as the classical trajectory does. So by using 
Eq. (|2~16l) one has 

\6xt(n > 0)| ~ - J== exp (?-e x n 3 / 2 ) . (2.25) 



The situation for the small modes is a bit complicated and requires careful considerations. As mentioned above the 
effective tachyonic mass of Sxk is of the order of e x H ^> H. As a result some modes become tachyonic even before 
horizon crossing. For modes which become tachyonic the spatial gradient term becomes negligible compared to the 
tachyonic mass and the evolution of <5%k becomes identical to that of the background solution. So for each mode it 
is important to find the time when it becomes tachyonic. Denoting the time when the mode k becomes tachyonic by 
n t (k), from Eq. (|2.2ip one has 

n t (fc)e 2 " t(fc) - f-^V • (2.26) 
\e x k c J 

Later on we shall call this time the "classicalization" time. The solution of the above algebraic equation are known 
to be given by the Lambert W function W(z), 

n t (k)= l -W(z)- z = 2^-^j , (2.27) 
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FIG. 1: Dynamics of the inflaton and waterfall fields (smoothed over a horizon-size patch). The solid blue line (the red 
squares) shows the numerical solution ( analytical estimate in Eq. (|2.10|l ) for inflaton field . The short dashed blue line 
for small N shows the dynamics of the inflaton field before the transition point. As one can see in this figure the inflaton 
dynamics before and after transition are connected smoothly. The dashed green line (the blue triangles) shows our numerical 
(analytical estimate) results for the dynamics of the waterfall field in a smoothing patch. In order to visualize this figure better 
the waterfall field amplitude is multiplied by a factor 2 x I0 3 . The model parameters are m = 6 x I0" 6 A'/p, s 2 = 3x 10" 8 , 
A = 0.14 and c/> c = M/g = 15M P , which give C = .15 and e x = 20. 



but we only need an approximate solution for nt(k) since nt(k) < 1 for the parameters of our interest. 
One can use the WKB solution Eq. (|2.23j) until n — nt(k), which gives 



<*x£(n) 



H 



2 k kp 



n < n t (k) . 



(2.28) 



After n = n t (k) this mode evolves as the classical background solution, Eq. (|2.16[) . By matching the solution with 
the classical trajectory at n = n t (k), for n > nj(fc) one approximately has 



*Xk(») 



H 



-e ™* exp 



; e x( n 



3/2 _ 3/2, 



n > n t (k) . 



(2.29) 



Here it is convenient to project Sxf i 71 ) & t time of the onset of waterfall n = as if all the modes were tachyonic at 
n > 0, as viewed in [29T |32|. This gives 



6x1(0) = 



H 



2/c ki 



exp 



-Wt(fc) - -e x n 3 t /2 (k) 



(2.30) 



Because of the exponential growth of the waterfall field quantum fluctuations after the transition, the e xpecta tion 
value (5x 2 ) becomes non-negligible and its rms value soon starts to behave as a classical field [2!| \ — \/ (&~X 2 h ■ As 



a result (^x 2 ) i s observed as a classical background for an observer within each Hubble horizon region [29j, |3lj, [3z 
Therefore it is necessary to calculate (5x 2 ) given in terms of its power spectrum Vs x as 



(tx 2 ) = 



d 3 k 
(2^)3 



d 3 k 
(2^)3 



dk 

T 



Vs x (k), 



(2.31) 



where the power spectrum is defined by 



(6XM = (2nrS%k + q)P x (k) , V Sx = ^P X (k) 



(2.32) 



7 



By using Eqs. (|2.25[) and (|2.29j) one can read off the power spectrum of the waterfall quantum fluctuations at n = 
as 



V Sx {k;0) = { 



H 2 ( k_ 



4tt 2 



H 2 e 2 



(2.33) 



The details of the analysis to calculate (8\ 2 ) is given in Appendix 1X1 where it is found that (Sx 2 ) is dominated by the 
small scale modes, 



(S X 2 (0)) ~ (8 X 2 (0)) t 



3e x /3 H 2 
16tt 2 



(2.34) 



As seen from Eq. (|2.33[) . there is a sharp peak in the spectrum of waterfall quantum fluctuations. In order to 
estimate the width of the peak, we expand the above spectrum around its peak. Solving dVs x /dn t = 0, the peak 
position is found as 



2e, 



2/3 



Expanding the spectrum around this momentum, one finds 



V Sx (k;0) 



4tt 2 



2/3 



with 



2ee, 



1 

3 \2Z 



exp 



2/3 



(n t (k) - n t (k max )Y 
2< 



(2.35) 



(2.36) 



(2.37) 



But we are interested in the width of spectrum in momentum space, <7*(k). By using Eq. (|2.26[) it is readily found as 

1 



cr*(/c) 



1 



2n t (k) 



\ + 0{e-^)\k n 



0.4fc„ 



(2.38) 
(2.39) 



This indicates that the width of the waterfall power spectrum is independent of the sharpness of the phase transition 
for large e x , which we verified also numerically. 



III. SN FORMALISM AND CURVATURE PERTURBATIONS 

In this section, using the 8 N formalism [Hcj, [H|, [l?} and [Hj], we calculate the curvature perturbations. In order 
to use the 6N formalism properly we trace back the number of e-folds from the end of inflation until the time of 
horizon crossing for each mode. To avoid confusion we denote the number of e-folds counted backward in time from 
the end of inflation by Af, that is, N = N e — N . Our strategy is to express N in terms of the fields <fi(n) and x 2 ( n ) 
(smoothed on every Hubble patch). 

For those modes which exit the horizon after the waterfall transition, by using Eq. (|2.11j) . one can easily find the 
curvature perturbation on comoving slices, TL C , as 

where 8<f> is to be evaluated on flat hypersurface as usual. Finding the curvature perturbation for those modes which 
exit the horizon before and during transition needs careful calculations. Here, we follow the same step as in [3lj for 
these modes. 
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Here it is worth mentioning that the duration of the waterfall stage is sensitive to the classical value of the waterfall 
field on every smoothing patch. As discussed before, because of the large tachyonic mass of the waterfall field there 
are modes whose spatial gradient can be neglected and which behave classically already before their scales cross the 
horizon, and hence affect the classical trajectory. Noting that the formula 1Z C = SAf is valid on scales over which 
small scale inhomogeneities can be smoothed out with a negligible influence on the geometry, we take the smoothing 
scale to be slightly larger than the comoving scale corresponding to the wavelength of the last mode which becomes 
tachyonic. 

As we discussed before, during the third inflationary stage (after completion of the phase transition) inflation 
proceeds in the form of chaotic inflation with a slight change in the effective mass of the inflaton. So the end of 
inflation is determined by the value of <f> alone as 



cf) = cf) e V2M P . 

From this point up to the time of the end of the waterfall transition, Af is given by Eq. (|2.f ip . 



Af = 



1 



AMp 



A" < A} , 



(3.2) 



(3.3) 



where Aff is the value of Af at the end of the waterfall transition. Here, for simplicity we consider that at the end of 
transition x is very close to its local instantaneous minimum so x 2 at A/" = Aff is 



X 2 (nf) = X 2 mm { n f) 



A£2 



(3.4) 



One can obtain Aff in terms of the number of e-folds from the critical epoch <j) — <f> c to the end of waterfall transition, 
rif, as 



A/>(n/) 



f 



4MJ, 



C 



by which one can readily find 

-W}>(i 

On the other hand using Eq. (12.91) one has 



2 & 



(3.5) 



(3.6) 



AMiSi 



■f 



By using the fact that 



^ 4> 2 C — AMptif one can find that 
5Af f (nf)~-5n f (l-2Cenf) 



dAff 



dn 



f 



-l + 2Cen/, 



(3.7) 



(3.8) 



in which e denotes the conventional slow-roll parameter which is approximately given by e ~ 2Mp/(f>^. 

Now we trace back the evolution to earlier times before the end of transition, Af > Aff. For this stage, instead of 
A/", it is more convenient to use n which is the number of e-folds from the critical point counted forward in time, i.e., 
n = rif +Aff — Af. Then x 2 ( n ) is given by 

x 2 (n) = exp [2 (f(n) - /(n/))] Xm in ( n f) > 
where f(n) for a sharp phase transition is given by 

f(n) = \e x . 



(3.9) 



(3.10) 



During this era, n is expressed in terms of (j>(n) as given by Eq. 



AMpn = <t>(n) 2 - $. 



1 - Cln -f- 



(3-11) 
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Here we note that n depends on n t and J\f in a non-trivial way, 

n(n f ,N) = n f +Aff(n f )-N'. (3.12) 
By virtue of the above geometric relation and by using Eq. (|3.8[) one has 

On 

— =2Cen f . (3.13) 
dn f 

Keeping in mind the above dependence of n on nj and TV, let us take the variation of Eqs. Q3.9p and (13. 1 1[) . Wc 
obtain 



5x 2 (n) 8Xmin( n f) 



(S X 2 (n)) 



-2M;5n 



+ 2f'(n)5n-2f'(n f )Sn f , 



C 



On the other hand from Eq. (|3.4[) . one finds 



(n)6(j>(n) ( 1 + — + Cen 



Xmin( n f) = AM P^T n f ■ 



A 

where we have used Eq. (12.91) and the fact that cf>f < </> c . This results in 

s xLin( n f) _ Sn f 



X 2 ■ 



i n f) 

Finally, solving Eqs. ((3~l"4l and (prT5]) for SAf, we find 

Sx 2 (n) dn 1 4>S(f>(n) 



5N = 



(S X 2 (n)) dn } -2r{n f ) + n-j l 



2Ml 



1 



C 
~2 



Ce 



2/'(n) 



~2f'(n f ) + nJ 1 dn f _ 



(3.14) 
(3.15) 

(3.16) 
(3.17) 

(3.18) 



To simplify the above relation we note that 2n$f(nf) 3> 1 which is valid in our model with a sharp phase transition. 
Then the above expression reduces to 



m = - 



Cerif Sx 2 {n) 
f'(n f ) (5x 2 {n)) 



l + j + Ce(n - 2n f ) 



2M 2 P 



(3.19) 



Noting that 8x 2 {n) and {x 2 { n )) have the same rt-dcpcndcnce the final result for the comoving curvature perturbation 
1Z C = SM can be expressed in terms of Sx 2 (0) and (x 2 (0))> 



n = 5Af = - 



S X 2 (0) 



f'(n f ) (6 X 2 (0)) 



1+ ( ± + C e(n-2n f ) 



2M\ 



(3.20) 



where and below we omit the suffix c from 1Z C and simply denote it by 1Z for notational simplicity. 

This is our key formula for computing the power spectrum and bispectrum. As can be seen from the above 
expression, the curvature perturbation has the conventional inflaton contribution up to slow-roll corrections in the 
inflaton mass and the contribution from the waterfall field. The latter is a dynamical effect which is intrinsic to our 
model, in contrast to many other models in which local features are added simply phenomenologically. 



IV. POWER SPECTRUM OF CURVATURE PERTURBATIONS 



Having found the final curvature perturbations with the SN formalism we now calculate the power spectrum. Our 
aim is to find an imprint of the sharp waterfall transition on the power spectrum. 

As clear from Eq. (|3.20j) . the power spectrum can be divided into two distinct contributions, the one from the 
inflaton and the other from the waterfall field, 



f' 2 (n f ) 



7>5 X 2 /X 2 + 



C 
2 



0{Ce) 



4Af4 



(4.1) 



Below we evaluate each contribution separately. 
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A. Contribution of the inflaton to the power spectrum 



Since the inflaton field is light throughout the whole stage of inflation the amplitude of its quantum fluctuations 
on flat hypersurface at the time of horizon crossing is given by the usual formula, 



6(/)(k) 



H(n k ) 



where H(nk) is the Hubble parameter at the time of horizon crossing. This gives 



4tt 2 



C 



I 2 ^H 2 



4Mf, 



(i + c) rv e + ff W 



48tt 2 



where V^(</>) is given by Eq. (|2.8|) . which may be rewritten to first order in C as 



Thus we obtain 



(i+cr (1- 



C <t> 2 \ m 



,2 jA 



(4.2) 



(4.3) 



(4.4) 



(4.5) 



Neglecting the corrections of 0(C), the above expression reduces to the standard formula, 'P^(fc) 
V 3 /(12TT 2 V' 2 M P )\ n=nk . 



B. Contribution of the waterfall to the power spectrum 



Now we calculate the contribution of the waterfall field perturbations to the power spectrum. In order to estimate 
this contribution we first note that 



((*X 2 ) k (*X 2 ) q ) = nv(fc)(2^ 3 (k + q). 



The correlation function of <5xk can be calculated using the identity [32| , 



m k (S X 2 ) q ) = 2 / ^\S X \^\W ^^(k + q) 



In Appendix [B] we show that 



V Sx >(h) 



^h 2 

~A^~ 2 



(4.6) 
(4.7) 

(4.8) 
(4.9) 



where £ is a numerical factor of order unity. This means that the power spectrum of Sx 2 is proportional to the power 
spectrum of 5\- Plugging this into Eqs. (|4.1j) and (I4.7[) and using the explicit form of f(n) given in Eq. (|3. 10[) and 
(<5x 2 (0)) calculated in Eq. (TAlTj) . we find 



fn ^_AC 2 e 2 n f eVs x (k;0) 



3 io/3 {Sx 2 {0)) 



(4.10) 



Finally, plugging the waterfall field perturbation spectrum (|2.33[) to the above, we obtain 

3 



16 r^2„ ,-20/3 > 3 ( k 

-Cen f e x £ I- 



16 



fC 2 e 2 n f e x 8/3 en t (k)e X p 



3/2 



3 c X"-t 



(k) 



k < c kf 



(4-11) 
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C. Total curvature perturbation spectrum 

We now consider the total curvature perturbation spectrum by adding contributions both from the inflaton field and 
the waterfall field. Since the waterfall contribution is peaked at k = k rnax > k c , let us first compare the amplitudes 
of V^(k) and V^(k) at k ~ k max . Using Eq. (|2.35j> for k max and comparing Eqs. (|4. 1 1|) and (|4.5|) . we find 

T> wf (k ) / f \ 4 ff \ -10/3 



10 J C 



where the approximation m/Mp ~ 10~ 6 has been used in order to satisfy the COBE normalization. This result shows 
that there can be a prominent peak even for a small C, say C ~ 0.1. 

In Fig. [5] we plot the total curvature perturbation power spectrum for the parameters C — 0.15, e = 0.01 and 
e x = 20. The peak at k — k max ~ fc c in the spectrum is due to the waterfall field contribution, (k). The spectrum 
away from the peak is dominated by the inflaton contribution, Vi(k). 

In passing, it is instructive to estimate nf, the duration of the phase transition. As mentioned before, we treat (<5x 2 ) 
as the averaged classical value of the waterfall field on each horizon size patch. A good criterion for the completion 
of the phase transition is when (^x 2 ) reaches the value of its local minimum Xmin( n f) gi yen by Eq. Q3.4p . With the 
help of Eq. (j3~TTj) . we find 

.2 A#2 



In Appendix IA1 the expectation value (Sx 2 (n)) is evaluated as 



( 5x \n f ))^-^n f . (4.13) 



(Sx 2 (n)) = (<5 X 2 (0)) exp(-e x n 3 / 2 ) ~ -^-^ exp(- £x n 3 / 2 ) . (4.14) 



^3 ' 16tt 2 " v 3 

Equating this with Xmin( n f)i we obtain an estimate 



n f 



r eZ 2/3 ; r ~ In 



x 



2 2/3 



327r^e x 
6A~ 



2/3 



(4.15) 



For our numerical example we find n/ ~ 0.5 so the phase transition is fairly sharp. But it is smooth enough to render 
the dynamics of the phase transition adiabatic. Namely, we are free from possible violations of the adiabaticity of the 
inflaton vacuum state that may occur for very sharp transitions as discussed in the literature HI, [l(| [ij], [l8|, . 

V. BISPECTRUM AND NON-GAUSSIANITIES 

Now we compute the bispectrum of this model. Due to the intrinsic non-Gaussian nature of Sx 2 , we expect to see 
large spiky non-Gaussianities when (k max ) > V^ikmax)- 

So far in our analysis, we have expanded SN up to Sx 2 as given by Eq. (|3.20[) . In order to calculate the bispectrum 
we need to expand SN up to Sx 4 - This is done in Appendix [C] The three point function can be read from Eq. (|C3I) as 

(ft kl ft k2 ft k3 ) = B 7 j(k 1 ,k 2 ,k 3 )(2^) 3 ( 5(k 1 +k 2 +k3) 

= (^) 3 <(^ 2 ) kl (<5x 2 ) k2 (^ 2 ) k3 ) 

+ \(N x ,fN x ^( [(Ax 2 )\ (Sx\ 2 (Sx\ 3 +cp) 

+ i(7V, ) 2 iV i00 ((^ 2 ) ki ^ k2 ^ k 3+c.p.), (5.1) 

where c.p. represents cyclic permutations, (ki,k2,k3) — > (k2,k3,ki) — > (k3,ki,k2) and A% 2 is the fluctuations of 
5x 2 (n,x) on scales larger than the horizon scale, as defined in Eq. (|2.14l) . 

There are two distinct contributions to the three point function. The first term in Eq. (|5.1I) is due to the intrinsic 
non-Gaussianity of Sx 2 ■ We define the intrinsic bispectrum of Sx 2 in the standard way by 

((Sx 2 )^X 2 )^X 2 ) k3 ) =i? 5x2 (fc 1 ,fc 2 ,fc3)(2^ 3 <5 3 (k 1 +k 2 +k 3 ). (5.2) 



12 



x 10 



, 9 



3 





-1 







2 



3 



4 



5 



lnk/k 



c 



FIG. 2: Power spectrum of the curvature perturbation. The dashed red curve shows an analytical estimate of the total 
curvature perturbation power spectrum by adding both the inflaton and waterfall field contributions, V^(k) + (k), given 
by Eqs. (|4.5|) and (|4.11[) . The blue solid curve shows the total curvature perturbation, Vn(k), in which the convolution integral 
Eq. (|4.8p is numerically calculated. As one can see in this figure, P™ / (k) peaks near k max ~ k c and decays quickly for k not 
close to kmax ■ The parameters are the same as in Fig. [1] 

The second term in Eq. (|5.1I) is due to nonlinear dynamics of the waterfall field, while the last term is that of the 
inflaton field which generically gives a negligible contribution when the inflaton is slow-rolling. Below we compute 
the first and second terms separately. 



Here we have intentionally avoided to call the above the dynamically generated bispectrum, because we shall see that 
it also includes some contribution from the intrinsic bispectrum of <5% 2 . 
By noting that 



A. Dynamically generated bispectrum 



Let us first concentrate on the second term, 



<^ kl ^ k2 ft k3 ) (2 ) = -{N^ fN^^ [(Ax 2 ) 2 ]^ (<5x 2 ) k2 (<5* 2 ) k3 +cp.). 



(5.3) 




(5.4) 



one has 




2 x 3 x (^ )ki {Sx%s ) + ( {SxXi (^2)^ ) + , p . 



6(S X 2 )B SX 2 (k u k 2 , k 3 ) (2ir) 3 S 3 (k 1 + k 2 + k 3 ) + ((5x 4 ) kl (^ 2 ) k2 (^ 2 ) k3 ) + c.p. 



(5.5) 



The first term on the right hand side above, proportional to B$ x 2, is of the same form as the first term in Eq. (|5.1I) . 
which we evaluate later. Here we first focus on the other terms in the square brackets. 
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We note that if 8\ 2 were Gaussian, we could express these terms in terms of the product of two point correlation 
functions. But since this is not the case for S\ 2 , there is also a contribution proportional to the three point correlation 
function of 5x 2 , as we shall see below. 

By expanding 5\ one has 

{(^ 4 ) kl (*X 2 ) k2 (S X \ 3 ) ((^ 2 ) kl _ q (*X 2 ) q (*X 2 ) k2 (*X 2 ) k3 ) , (5.6) 



in which we have introduced d 3 q — d 3 q/(2ir) 3 to simplify the notation. To calculate the r.h.s of this equation we 
should first classify possible contractions. Since we are not interested in tad-pole type graphs but only in irreducible 
graphs, there are some restrictions on non-trivial contractions. First, contractions between the terms within any 
of (<5x 2 )p themselves are not allowed. Second, contractions should not be closed only within the terms in (<5x 4 ) kl , 
corresponding to the first two terms in the r.h.s. of the equation. Then one finds there are two different classes of 
contractions. The first class is in which there is one contraction between a pair of terms in (<5x 4 ) kl . The second class 
is in which there is no contraction between any pair of terms in (#x 4 ) kl themselves. 

Let us consider the first class and count the number of possible contractions. There are 4 choices to choose a pair 
in (5x 4 )ki- Then one of the remaining two S\ has 4 choices to contract with one of 8\ in (<$X 2 )k 2 an d (^X 2 )k 3 , and 
the last 5x in (<5x 4 )ki has 2 choices to contract the remaining terms in (^x 2 )k 2 an d (^X 2 )k 3 - Thus there are in total 
4 x 4 x 2 = 32 possible contractions in the first class. They all give the same result. So let us calculate one of them: 

F===\ 



J d 3 q~[[ d 3 Pl / (^Xki-pi-q^Xpi) (^Xq-P2^Xp 2 ) ( <5 Xk 2 -p 3 ^Xp 3 ) (<5Xk 3 - P4 <5X] 

= J'd J qY[d 3 pi\Sx\ cl -p 2 \\ 2 \Sxp 1 \ 2 \SXp2\ 2 \SXp a \ 2 * 5-factor, (5.7) 

i 

where 

^-factor = <5 3 (ki - Pl - p 2 )<5 3 (k 2 + pi - p 3 )£ 3 (k 3 + p 2 p 4 )<5 3 (p 3 + P4) ■ (5.8) 
Performing first the integrals over p 2 , p 3 and P4, and then over q, the above reduces to 

J d 3 q J d 3 Pl |5x| q+Pl ^k 1 || 2 |5x Pl | 2 |<5x|k 2 +p 1 || 2 |'5x|k 1 -p 1 || 2 '5 3 (k 1 +k 2 +k 3 ) 

= (5 X 2 ) J d 3 Pi|5x Pl | 2 |5x|k 2+ p 1 || 2 |'5x|k 1 -p 1 || 2 ( 5 3 (k 1 +k 2 +k 3 ) 

= J^ 2 )((^ 2 ) kl (^ 2 ) k2 (^ 2 ) k3 > ) (5-9) 



where we have used the fact that 3z 



((^M^X^X^ 8 / ^IJxpJ'l^lfe+pxiH^lkx-pxll^^x+ka + ka). (5.10) 



Thus we obtain 



J d 3 qY[d 3 p l ((Sx\, 1 -p 1 -qSxp 1 )(Sxq-p 2 Sxp 2 ) (^Xk 2 -p3^Xp 3 ) (foks-p^Xp*)) 

i 

= i(5x 2 )S 5x2 (fc 1 ,fc 2 ,fc 3 )(27r) 3 <5 3 (k 1 +k 2 + k 3 ), (5.11) 

where the bispectrum Bg x 2 {k\, fc 2 , fc 3 ) is defined in Eq. (I5.2p . Since there are 32 x 3 of the same terms, the contribution 
from this class amounts to 



<(^X 4 )k 1 (5x 2 )k 2 (5x 2 )k 3 )+c.p. 



1 

* ~' '''S X 

1st 



32 x 3 x -(Sx 2 }B gx 2(k u k 2 , k 3 ) (2 7 r) 3 J 3 (k 1 + k 2 + k 3 ) 



12 (Sx 2 )B Sx 2(ki,k 2 , k 3 ) (2^) 3 «5 3 (k 1 + k 2 + k 3 ) . (5.12) 
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Let us now turn to the second class, in which every Sx in (#x 4 )ki is contracted with one of <5% i n either (^x 2 )k 2 or 
((5x 2 )k 3 - Hence, there are 4 x 3 x 2 = 24 possible contractions. Again they all give the same result. One of them is 
given by 



d 3 <7jJd 3 K / (gXki- Pl -q^X Pl ) (SXq-pJXpi) (^Xk 2 -p 3 ^X P3 ) (foka-p^XpJ 



d 3 qY[d 3 pi \5x\k 2 -p 3 \\ 2 \Sx Pl \ 2 \$Xp2\ 2 \dXp 3 \ 2 x 5-factor, 



where 



5-factor = S 3 (k 1 + k 2 pi p 3 q)<5 3 (k 3 + pi - P4)<5 3 (q - P2 + P 3 )<5 3 (p2 + P4) • 
Performing first the integrals over p 2 , P3 and P4, and then over q gives 



(5.13) 



(5.14) 



d 3 Pi \s Xpi ri*xik»+ P1 ir / d d 9|5xiq+p 1 -k 1 iri^xiq+p 1 +k 3 ir 5 ( k i + ^ + k s) 

= ip, 5x2 (fc 2 )P (5x 2(A :3 )(27r) 3 ^(k 1 +k 2 + k3), (5.15) 

where we have used the definition of Ps x 2 , Eq. (|4.6I) . Since there are 24 of them, plus cyclic permutations, the 
contribution from this class is in total, 



<(^ 4 )k 1 (<5x 2 )k 2 (5x 2 )k3)+cp. 
1 



= 24 x 



= 6 



J 2nd 

[P SX 2 (k 2 ) P Sx 2 (fc 3 ) + c.p.] (2 7 r) 3 5 3 (k 1 + k 2 + kg) 
^(fci)^(fc2) + cp.] (2 7 r) 3 ,5 3 (k 1 +k 2 + k 3 ). (5.16) 
Adding up all the contributions given by Eqs. (|5.12j) and (|5.16[) together with the first term in (|5.5j) . we obtain 
([(A x 2 ) 2 ] ki (^ 2 ) k2 (<5 x 2 ) k 3+c.p.) 

= Q({5 X 2 )B Sx 2(k 1 ,k 2 ,k 3 )+ [iV(*i)iV(*2)+c.p.]) (2^) 3 «5 3 (k 1 +k 2 + k 3 ). (5.17) 

Substituting the above in Eq. (|5.3I) . the three point correlation function of the curvature perturbation can be repre- 
sented in terms of the power spectrum and bispectrum of Sx 2 as 



(^ki^k 2 '^k3)(2) — 



3(iV, x2 ) 2 iV ix2x2 ((5x 2 >^(fci,A :2 ,fc 3 ) 



+3(A^ x 2) 2 7V x2x 2 (p s%2 P^ (ha) + cp.) 



27r) J <5 3 (ki +k 2 + k 3 ) 



(5.18) 



As it is clear in the above result, while the second term is an ordinary nonlinear interaction term with the vertex 
proportional to N X 2 X 2 , which can be easily evaluated, the first term is due to the intrinsic non-Gaussianity of the Sx 2 
field which needs to be computed. 

The above equation can be further simplified by using Eq. ([C4| and noting that P% f = N 2 2 P Sx 2 , 



(TC kl ft k2 ft k3 ) (2) = -3(A^ iX 2) 3 P 5x2 (/c 1 , k 2 , fc 3 )(27r) 3 <5 3 (k x + k 2 + k 3 ) + (TZ^Tl^TZ^) dyn , 
where we have defined the dynamically generated bispectrum of the curvature perturbation by 
(TZ kl TZ^n k3 ) dyn = P^"(k 1 ,k 2 ,k 3 )(2 7 r) 3 ( 5 3 (k 1 +k 2 +k 3 ) 



(5.19) 



, N ,x 2 x 2 



(p- / (fc 1 )P- / (fc 2 )+c.p.)(27r) 3 <5 3 ( 



ki 



k 3 ) 



(5.20) 
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B. Bispectrum from intrinsic non-Gaussianity 

Now we evaluate the bispectrum from the intrinsic non-Gaussianity of <5% 2 . From Eqs. (|5.ip and f|5 . 19[) . we have 
(K kl K k2 1l k3 } in t = B™ t (k 1 ,k 2 ,k 3 )(2^) 3 ( 5(k 1 +k 2 +k3) 

= -2(^ x 2) 3 B 5x 2(ki,k 2 ,k 3 )(2^) 3 ( 5 3 (ki + k 2 +k3) . (5.21) 
Following [32^ . we obtain an expression for the bispectrum of Sy 2 as 

B Sx 2(k 1 ,k 2 ,k 3 ) = 8 /d^l^l^Xl^-qil^Xlk.+qil 2 . (5.22) 



It is hard to calculate the above integral in general, but it may be evaluated in the squeezed limit, say k 3 <C fci = k 2 = k. 
In this limit, the above reduces to 



B% 3 (k) = B Sx '(k,k,0)~8 / Wx 8 | 2 |*X|k- q |l 



In Appendix [B] this integral is evaluated, and the result (given by Eq. (|B14[) 1 is 



t/3 H 4 



1 



p- P <5x 2 ( fc ) 



(5.23) 



(5.24) 



where the second step follows from Eq. (|B10I) or (I4.9P , and £ and £' are constants of order unity. 

Another limiting case of interest is when the magnitudes of all the momenta are equal to each other, k\ = k 2 = k 3 , 
the so-called equilateral limit. Although we have no clue whatsoever to evaluate the bispectrum in this limit, it is 
plausible that the amplitude is at most of the same order as that in the squeezed limit, if not much smaller. So let us 
set 



where £ cq is a non-dimensional factor supposedly of order unity. 



(5.25) 



C. Total bispectrum and }nl parameter 



It is customary to express the bispectrum of the curvature perturbation in terms of the non-Gaussianity parameter 
Jnl [53]. The standard definition of the non-Gaussianity parameter in Fourier space is [Hf 



6, , U s B K (k 1 ,k 2 ,k 3 ) 

■7JNL{ki,k 2 ,k 3 ) = . 

5 [Pn{k 1 )P n {k 2 ) + c.p. 



In our case we may decompose it as 



TNL = Jnl 



fint _,_ fdyn , f(j> 

Jnl ' Jnl ' J, 



NL ' 



where 



^ pint 
~^JNL 



rdyn 
JNL 



Jnl 



B^{k ly k 2 ,k 3 ) 
[Pn(ki)Pn(k 2 ) + c.p. 

B^ n (ki, k 2 , k 3 ) 
[Pn{kx)Pn{k 2 ) + c.p. 



= -2 (N 



= 3 



, B Sx 2(k 1 ,k 2 ,k 3 ) 
[Pn(ki)Pn{k 2 ) + c.p.] 

F^ / (fc 1 )P^ / (fc 2 )+c.p. 



Pn(ki)Pn(k 2 ) + c. v . 



B^{ki,k 2 ,k 3 ) 


1 JV^ 


F^(fc 1 )F^(fc 2 )+c.p. 


[Pn(h)Pn(k 2 ) + c.p.] 


2 (^) 2 


i^(*i)^(*2) + c.p. 



(5.26) 
(5.27) 

(5.28) 
(5.29) 

(5.30) 
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As noted before, the inflaton contribution ffj L is known to be at most of the order of the slow-roll parameters, 
fffL = ^C 6 ' 7 ?) 0] hence can be safely neglected. So we focus on the contributions from the waterfall fields. 

First let us consider /^™. As clear from its form, it is non- negligible only when the amplitude of P^ (k) is 
comparable to or greater than that of P^(fc), and this may happen only at and around the peak of the spectrum 
k = k max . Then it is easy to see that fj^£ is non- negligible only when all of k\, ki and k-$ are approximately equal 
to k max . Thus setting k\ = k^ = k% = k and assuming P^ (k) dominates the spectrum, we find 



-f dyn (k-k ) ~ 3 'X 2 * 2 - 3e * ff5 sFi 

5 ( N ,x 2 ) Cen f 



where we have used Eq. (|C4I) in Appendix [Cl to eliminate N iX 2 and N X 2 X 2 from the intermediate expression. With 
e x 3> 1, C 1 and rif ~ 1, we see that /|%™ can become very large, centered around fc = k max . For example, for the 
parameters used in our numerical analysis, we have 



Next we consider As clear from the form of B,5 X 2 in Eq. (I5.21j) . or its explicit evaluation, it remains finite in 

the squeezed limit, while the denominator of f™[ diverges since P^ is approximately proportional to k~ 3 . Therefore, 
fj^l is completely negligible in the squeezed limit. In contrast, the denominator is finite and of order k~ 6 in the 
equilateral limit. Hence, using Eq. (15.25[) , we obtain an estimate 

5 f NL (k)-S 3fc 3 p 2 (fc) N tX , 3 klv ^ {k) ■ (5-33) 

Again, this contribution to the non-Gaussianity is exponentially negligible except at around the peak of the waterfall 
field spectrum. Hence assuming P^J ~ "P^, and manipulating with the help of Eqs. (|4.10l) . f|4. and (|C4I) . the above 
estimate gives 

int c H 2 ^ 3 1 2 dyn 

fNL(kmax) ~ ^° q ^ x 2 £3 P n (k) ~ ^ C<le ^^ NL ^ max ^ ' (5.34) 

Depending on the value of the parameter ^ cq of which we have no quantitative estimate, the intrinsic non-Gaussianity 
can be larger than the dynamical non-Gaussianity. 

In any case, we conclude that the total non-Gaussianity parameter Jml is at least as big as 



fNL(kmax) — PnL + Inl ~ r< * ' (5.35) 



The width of this sharp feature in the bispectrum is the same as that of the spectrum, Eq. (I2.38|) . ie, a{k) ~ 0Ak ri 



VI. CONCLUSION AND DISCUSSIONS 



In this paper we presented a dynamical mechanism to generate a sharp feature during inflation. The key role is 
played by the waterfall field x which becomes tachyonic during inflation. We work in a region of the parameter space 
where the waterfall transition is fairly sharp, i.e., the duration of the transition is about one e-fold or so. Because of 
the coupling g 2 <f> 2 x 2 the phase transition induces a sharp but small change in the inflaton mass. In much of previous 
works in which sharp changes in the inflaton mass were studied, the focus was on the dynamics of the inflaton itself. 
In contrast, in our model a local feature is induced by non-trivial interactions between the waterfall and inflaton fields. 
In particular, the waterfall quantum fluctuations played the key role in determining the local feature. 

Before the phase transition x is very massive, so it has no classical evolution. It stays at its local minimum at x — 0. 
When the inflation passes a critical value, x becomes tachyonic and the waterfall transition commences. Then the 
squared fluctuation, (5x 2 ), starts to grow exponentially and the effective classical trajectory is determined by (Sx 2 ) 
averaged over each Hubble horizon patch. The fluctuations from one horizon region to another, Ax 2 = Sx 2 — (^X 2 )j 
determines the fluctuation around the classical trajectory. 

We calculated the power spectrum of A% 2 which is found to have a peak near the comoving scale k c that crosses the 
horizon at the onset of transition, k = k max ~ k c . This in turn induces the curvature perturbation which shows a peak 
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near k = k ma x- The ratio of thus induced curvature perturbation to the conventional curvature perturbation from 
the inflaton field depends on the model parameters as given in Eq. (I4.12[) . For reasonable values of the parameters, 
we find this ratio can become of order unity or even larger than unity. 

This local feature we found may be used to explain the glitches found in the observed CMB angular power spectrum. 
One may also consider many waterfall fields coupled to the inflaton to produce a series of waterfall phase transitions 
and induce multiple local features in the curvature perturbation. It would be interesting to study numerically the 
effect of single or multiple local features in our model and compare it with CMB data. 

Due to intrinsic non-Gaussian nature of 5x 2 distribution one expects to see large spiky non-Gaussianities [551 ] when 
(k max ) becomes comparable to V^. We have shown that Jnl has both intrinsic and dynamical contributions. The 
intrinsic contributions originates from <5% 2 three-point function while the dynamical part comes from the non-linear 
dynamics of the waterfall field. It is shown that fNL{k m ax) ~ e x/ e ^' m wm ch e x measures the sharpness of the 
waterfall phase transition. As a result, the sharper is the phase transition the larger is Jnl- As in the case of power 
spectrum, the bispectrum is highly peaked at k ~ k max . It would be very interesting to investigate the observational 
consequences of these spiky non-Gaussianities. 
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Appendix A: Variance of waterfall field quantum fluctuations 

Here we calculate (6x 2 (n)) in some details. As discussed in the text, we must consider this as part of the classical 
background after it begins to evolve as a classical field, 



exp 



3 x 



(Al) 



Then it is convenient to define (8x 2 (0)) not by its actual value at n = 0, but the value it would take if it evolved 
classically from the beginning. That is, for modes that become tachyonic by the end of the waterfall transition n — n/, 
we define 



<<5 X 2 (0)) = (S X 2 (n f ))cxp 



4 3/2 

3 x / 



(A2) 



As mentioned in the text, we divide it into small scale and large scale parts, denoted by subscript S and L, respectively, 



(6x 2 m = (6x 2 (0))s + (6 X 2 (0))i 



(A3) 



As shown in [321 ] . for the large scale part we have 



(<5x 2 (0))l = 



H- 



(A4) 



As for the contribution of small scale modes, as we show below, it was somewhat over-estimated in the previous 
literature. A more accurate result is obtained as follows. 

As discussed around Eq. (|2.29[) in the text, we approximately have 



SxUn) 



H 



■e "* exp 



1 I 3/2 3/2^ 

-e x (n ' -n t < ) 



n > n t (k) . 



(A5) 



Literally speaking, however, the above does not reproduce the correct behavior when n — n t -C 1. At this stage, it 
behaves as 



H 



■e "* exp 



\t x {n-n t fl 2 



(A6) 
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If we simply extrapolate this to n = n/, we obtain 

H 



e ™* exp 



-6^/ - 7l t ) 3 / 2 



(A7) 



instead of the one that follows from the approximate formula (|A5 



H 



■e ™* exp 



2 , 3/2 3/2^ 

-e x {n/ -n t ' ) 



(A8) 



Then the formula (IA7|) gives an estimate of $xf (0) as 



6x1(0) = 



H 



,5xf(0)exp 



e n * exp 



-e x (n / - n t ) 3/2 ~ e x n 



3/2 



e x (n/ - n t ) 3/2 - e x (n 3 / 2 - nj 2 ) 



(A9) 



where Jxf (0) is estimated by using Eq. (|A5[) . We easily see that Sxi (0) < ^xf (0)- Thus one would expect that using 
$Xk (0) would give a slight underestimate, while using <5xf (0) would give a slight overestimate. 
Let us first estimate (<5x 2 (0))s by using <5xf (0). Using Eq. (|2.26[) . we find 



2 fr2 



(^ 2 (0))5 = 



4ir 2 



«t=0 



dn t (n t + -)e 



(A10) 



— 2/3 

The exponential dependence of the integrand introduces a natural cut-off at n cut ~ e x . Hence one can neglect nt 
with respect to 1/2 in the integrand, and extend the integral range to infinity to obtain 



e 2 x H 2 



8tt 2 



elH 2 



4 /3 h 2 



dn* e -4 / 3 ^™' = — 75 



(All) 



This gives a slightly overestimate of the true (5x 2 (0))s- 

Now let us consider using <5xf (0). Instead of the integral (|A10[) . we have 



n ( 2 ,H 2 
(6x 2 (0))s = ^ r 



1 1 
dn t (n t + -) exp 

n f =0 * 



4 3/2 



i - 



77, 



3/2 



(A12) 



Again, similar to the previous case, the integral is dominated by the integrand at nt/nj <C 1. Then expanding the 
exponent in n t /n/, we may approximate it by 



<<5x 2 (0)K< 



e 2 H 2 r n f 



8tt 2 



dnt exp 



n t =0 



1 1/2 



e 2 x H 2 I 



8?r 2 2e x n/ 2 Wtt 2 ,,] 



1/2 
/ 



(Al3) 



This result is by a factor e^n^ 2 smaller than the estimate (IAll[) . For a typical value of e x and nj, say e x ~ 10 

^ /3 l /2 

and n/ ~ 0.5, however, this factor is e x ~ 1.5. So qualitatively there is not much difference between the two 

estimates. 

To summarize, we conclude that the dominant contribution to the variance of the small scale waterfall field quantum 
fluctuations comes from the modes around e^n^ 2 < n < e x 2 ^ 3 , and given by somewhere between Eqs. (| Al 1|) and 
(|A13[) . Comparing these with Eq. (|A4[) . we readily see that the dominant contribution to the variance of the waterfall 
field quantum fluctuations comes from the small scale modes, 

(Sx 2 (0)) = (Sx 2 (0))s + (S X 2 (0))l * (Sx 2 (0))s ■ (A14) 
Since there is not much difference between Eqs. (jAllI) and (IA13[) . for definiteness we use Eq. (|A11[) in the text. 
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Appendix B: Correlation functions of square of waterfall field quanta V8\ 2 

This appendix is devoted to find a good approximation for the correlation functions of the 8\ 2 appearing in the 
power spectrum and bispectrum analysis, Eqs. (|4.8[) and (|5.22[) . 

First, we work out the following convolution integral which is necessary is calculations of power spectrum of the 
waterfall field 

(MkM 9 ) = 2 / (|^N k - q || 2 |^| 2 (2-) 3 <5 3 (k + q). (Bl) 
The above integral can be divided into the radial and angular parts as follows 

2 J -0^\8x\^\\ 2 \S Xq \ 2 = J dn q V Sx (q) J d(- cos0)|5 X |k- q | | 2 . (B2) 

The above integral can not be calculated analytically but can be estimated as follows. First, note that as the 
power spectrum of waterfall field Vs x is highly peaked around k = k max so one can expect that the peak of curvature 
perturbation V-r, ~ V$ x -2 also take place near k max 3> k c . This assumption will be checked in the following. The other 

main point is that Eq. (|2.33[) shows that |#X|k-q|| i s constant value wA for |k — q| < k c and exponentially decay for 
|k — q| > k c . As a result one can conclude that the second integral in Eq. (|B2[) is negligible except for 

|k - q| < Hk c (B3) 

in which the numerical factor £ can be estimated from Eq. (|2.30[) . To find an estimate of £ it is natural to estimate 
the width of integration, Afc, the point at which P$ x decreased by a factor 1/e from its value at k = k c , Ps x (k c ) = 
i?Q/(2fc 3 ). So by using Eq. (I2.30|) and noting that for k ~ k c , n t (k) <C n*(/c), one has 

£ ~ e. (B4) 

In order to satisfy condition (|B3[) one can simply find that the amplitude of integral momentum q should be near k 
and at the same time the angle between k and q, denoted by 9, should be near zero. By defining 

q = k + Aq, Aq < k, and 6 = + AO (B5) 

one can find that the condition 

|k - q| 2 < ek 2 c (B6) 

results in 

Ag<£fc c — > An ^^Y (B7) 
A(-cosO) = (B8) 
With these approximations the integral Eq. (|B2[) can be estimated for k 3> k c as 

2 J (03l^|k- q || 2 |^ 9 | 2 = 2 I dn q d(~ cos 9)\S X \ k -^\ 2 V lSxl (q) 

~ An q A(- cos 6) ^Vs x (k) 

= ^^x(fc) (B9) 
Finally one can find the following good approximate for Vs x 2 in terms of Vs x 



P tfW-^^xW' (BIO) 



Similarly, we present an approximation for the three point correlation function of 5\ 2 in the squeezed form 

((^ 2 ) kl ('5x 2 ) k2 (^ 2 ) k3 ) sq ^8 ( 5 3 (k 1 +k 2 +k3) J d\\5 Xq \ 2 \5 X ^. (Bll) 
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One can approximate the above integral in the same fashion as we did above for the two point correlation function 

((<5 X 2 ) kl (S X 2 ) k2 (<5x 2 ) k3 ) sq ^8(27r) 3 < 5 3 (k 1 +k 2 +k 3 ) J dn q V Sx (q) J d{- cosfl)|<5 X |k- q | | 4 (B12) 
By using the above relation, one finds 

dn q Vs x (q) J d(- cos 0)|5x|k-q| I* 

El- 
V 6 ' 



R sq 



8An q A(~co S 0)-^V Sx (k) 



k 3 k 3 



V Sx {k), 



in which the numerical factor similarly defined as £, is obtained to be £' ~ y^e. Finally one has 

^(*)*f£||Ux(*)- 



(B13) 



(B14) 



Appendix C: 8M up to order A^ 4 

In order to find 8J\T up to order A^ 4 it is enough to extend Eq. (|3.14[) to next order in Sx 2 , which reads 
A X 2 (n) 1 A X 4 (n) _ S X 2 mm (n f ) , 9f ,, Un 



This modifies SAf in Eq. (|3~18| to 

(57V = 



A X 2 (n) 1 A X 4 (n) 



(S X 2 (n)) 2 ( X 2(n))2 



4> 8<f>(ri) 
2M 2 p 



C 

1 + — + Cen 



For a sharp transition, rif < 1, the above expression reduces to 



«/V = - 



Cen/ / A X 2 (n) 1 A X 4 (r 



f'(n f ) \(S X 2 (n)) 2<rW 



9n / -2/'(n / )+nJ 1 

2/'(n) 9n 
^/'(n^ + n/dn/ 



1 + - + Ce(n - 2n/) 



2MJ, 



From the above, we find the relation between -/V x 2 and N X 2 X 2 as 



n x , = -(s x 2 (o))n x . 



-Ct 



f'(n f ) (6 X 2 (0)) ~ n 1 / 2 e x (6x 2 (0)) 



(CI) 



(C2) 



(C3) 



(C4) 



where the second equality follows from the fact that f'{n) = e x n 1 / 2 . 
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